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ABSTRACT 

The purpose of this paper is to explore the dynamical behaviour of hot accretion flow 
with thermal conduction. The importance of thermal conduction on hot accretion flow 
is confirmed by observations of the hot gas that surrounds Sgr A* and a few other 
nearby galactic nuclei. In this research, the effect of thermal conduction is studied 
by a saturated form of it, as is appropriate for weakly collisional systems. The an- 
gular momentum transport is assumed to be a result of viscous turbulence and the 
a-prescription is used for the kinematic coefficient of viscosity. The equations of accre- 
tion flow are solved in a simplified one-dimensional model that neglects the latitudinal 
dependence of the flow. To solve the integrated equations that govern the dynamical 
behaviour of the accretion flow, we have used an unsteady self-similar solution. The 
solution provides some insights into the dynamics of quasi-spherical accretion flow and 
avoids from limits of the steady self-similar solution. In comparison to accretion flows 
without thermal conduction, the disc generally becomes cooler and denser. These prop- 
erties are qualitatively consistent with performed simulations in hot accretion flows. 
Moreover, the angular velocity increases with the magnitude of conduction, while the 
radial infall velocity decreases. The mass accretion rate onto the central object is re- 
duced in the presence of thermal conduction. We found that the viscosity and thermal 
conduction have the opposite effects on the physical variables. Furthermore, the flow 
represents a transonic point that moves inward with the magnitude of conduction or 
viscosity. 
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' 1 INTRODUCTION 

\ Accretion is one of the most important physical processes in 
astrophysics. It is widely accepted that the accreting matter 
toward the central object is a source of power active galactic 
nuclei (AGN), and galactic X-ray sources (see for a review 
e.g. Frank et al. 2002). This idea is also well-applicable to in- 
terpret many observations of astrophysical phenomena, such 
as, prototype stellar objects (Chang & Choi 2002), symbi- 
otic stars (Lee & Park 1999), gamma-ray bursts (Brown et 
al. 2000). There are some types of accretion flows that a 
significant fraction of the generated heat by dissipation pro- 
cesses retains in the fiuid rather than being radiated away, 
have been the subject of considerable attention in recent 
years (Ogilvie 1999; Chang 2005, Akizuki & Fukue 2006, 
KhesaU & Faghei 2009, Faghei 2011). These advection dom- 
inated accretion flows (ADAF) place an intermediate po- 
sition between the spherically symmetric accretion flow of 
non-rotating fiuid (Bondi 1952) and the cool thin disc of clas- 
sical accretion disc theory (e. g. Pringle 1981). These types 
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of accretion flows have been widely applied to explain obser- 
vations of the galactic black hole candidate (e. g. Narayan 
et al. 1996, Hameury et al. 1997), the spectral transition of 
Cyg x-1 (Esin 1996) and multi- wavelength spectral proper- 
ties of Sgr A* (Narayan & Yi 1995; Manmoto et al. 2000; 
Narayan et al. 1997). 

The X-ray observations of black holes imply that they 
are capable of accreting gas under a variety of flow config- 
urations. In particular, observational evidences confirm ex- 
istence of hot accretion flow, contrasted with the classical 
cold and thin accretion disc scenario (Shakura & Sunyaev 
1973). Hot accretion flow can be found in the population of 
supermassive black holes in galactic nuclei and during qui- 
escent of accretion onto stellar-mass black holes in X-ray 
transients (e.g., Narayan et al. 1998a; Lasota et al. 1996; 
Di Matteo et al. 2000; Esin et al. 1997, 2001; Menou et 
al. 1999; see Narayan et al. 1998b; Meha & Falcke 2001; 
Narayan 2002; Narayan & Quataert 2005 for reviews). Chan- 
dra observations provide tight constraints on the density 
and temperature of gas at or near the Bondi capture radius 
in Sgr A* and several other nearby galactic nuclei. Tanaka 
& Menou (2006) used these constraints (Loewenstein et al. 
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2001; Baganoff et al. 2003; Di Matteo et al. 2003; Ho et al. 
2003) to calculate the mean free path for the observed gas. 
They suggested that accretion in these systems will be pro- 
ceeded under the weakly collisional condition. Furthermore, 
they suggested that thermal conduction can be as a possible 
mechanism by which the sufficient extra heating is provided 
in hot advection dominated accretion flows. 

Generally, semi-analytical studies of hot accretion flows 
with thermal conduction have been related to steady state 
models (c. g. Tanaka & Mcnou 2006; Johnson & Quataert 
2007; Shadmehri 2008; Abbassi et al. 2008, 2010; Ghanbari 
et al. 2009), and dynamics of such systems have been stud- 
ied in simulation models (c. g. Sharma ct al. 2008; Wu ct al. 
2010). For example, Tanaka & Mcnou (2006) have carried 
out a related analysis and found the accretion flow can spon- 
taneously produce thermal outflows driven in part by con- 
duction. Their analysis is two-dimensional but self-similar in 
radius. Their assumption of self-similarity enforces a density 
profile that varies as r^^^'^, whereas simulations of ADAFs 
consistently find density profiles shallower than this (e.g., 
Stone et al. 1999; Igumenshchev & Abramowicz 1999; Stone 
& Pringle 2001; Hawley & Balbus 2002; Igumenshchev et al. 
2003). Johnson & Quataert (2007) studied the effects of elec- 
tron thermal conduction on the properties of hot accretion 
flows, under the assumption of spherical symmetry. Since, 
electron heat conduction is important for low accretion rate 
systems, thus their model is applicable for Sgr A* in the 
Galactic centre. They show that heat conduction leads to 
supervirial temperatures, implying that conduction signif- 
icantly modifles the structure of the accretion flow. Their 
model similar to Tanaka & Menou (2006) was the steady 
state, but they solved their equations numerically. 

As mentioned, semi-analytical studies of hot accretion 
flows with thermal conduction have been in a steady state. 
Thus, it will be interesting to study dynamics of such sys- 
tems. Ogilvie (1999) by the unstcadj' self-similar method 
studied time-dependence of quasi-spherical accretion flow 
without thermal conduction. The solutions of Ogilvie (1999) 
provided some insight into the dynamics of quasi-spherical 
accretion and avoided many of the limits of the steady self- 
similar solution. In this research, we want to explore how 
thermal conduction can affect the dynamics of a rotating 
and accreting viscous gas. We answer this question by solv- 
ing Ogilvie (1999) model that is affected by thermal con- 
duction. This paper is organized as follows. In Section 2, we 
define the general problem of constructing a model for hot 
quasi-spherical accretion flow. In Section 3, we use the un- 
steady self-similar method to solve the integrated equations 
that govern the dynamical behaviour of the accreting gas, 
and numerical study of the model is brought in this section, 
too. We will present a summary of the model in Section 4. 



2 BASIC EQUATIONS 

We start with the approach adopted by Ogilvie (1999), who 
studied quasi-spherical accretion flows without thermal con- 
duction. Thus, we derive the basic equations that describe 
the physics of accretion flow with thermal conduction. We 
use the spherical coordinates (r, 9, <p) centred on the accret- 
ing object and make the following standard assumptions: 



terized by the Newtonian potential of a point mass, *t = 
— GMt/r, with G representing the gravitational constant 
and M« standing for the mass of the central star. 

(ii) The written equations in spherical coordinates are 
considered in the equatorial plane 9 — 71/2 and terms with 
any 9 and 4> dependence are neglected, hence all quantities 
will be expressed in terms of spherical radius r and time t. 

(iii) For simplicity, self-gravity and general relativistic ef- 
fects have been neglected. 

Under these assumptions, the dynamics of accretion 
flow describes by the following equations: 
the continuity equation 



dp 1 i9 / 2 

the radial force equation 
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the azimuthal force equation 
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Here p the density, Vr the radial velocity, Q the angular ve- 
locity, 0,k[= {GMa,/r^y^^] Keplerian angular velocity, p the 
gas pressure, 7 is the adiabatic index, u the kinematic vis- 
cosity coefficient and it is given as in Narayan & Yi (1995a) 
by an a-model 



Pgas 
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The parameter of a is assumed to be a constant less than 
unity. The terms on the right-hand side of the energy equa- 
tion, Qvis is the heating rate of the gas by the viscous dis- 
sipation, Qrad represents the energy loss through radiative 
cooling, and Qcond is the transported energy by thermal con- 
duction. For the right-hand side of the energy equation, we 
can write 



Qa 



Qv 



Qrad H" Qcond 



(6) 



where Qadv is the advective transport of energy. We employ 
the advection factor, / = 1 — Qrad/Qvis, that describes the 
fraction of the dissipation energy which is stored in the ac- 
cretion flow and advcctcd into the central object rather than 
being radiated away. The advection factor of / in general 
depends on the details of the heating and radiative cooling 
mechanism and will vary with position (e.g. Watari 2006, 
2007; Sinha et al. 2009). However, we assume a constant / 
for simplicity. Clearly, the case f = 1 corresponds to the 
extreme limit of no radiative cooling and in the limit of ef- 
flcient radiative cooling, we have / = 0. 

In a collisional plasma, mean free path for electron en- 
ergy exchange. A, is shorter than temperature scale height, 
Lt = T/|VT|, and thus the heat flux due to thermal con- 
duction can be written as 



(i) The gravitational force on a fluid element is charac- 



Fcond = -kVT, 



(7) 
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where k is the thermal conductivity coefficient. Thermal 
conductivity in a dense, fully ionized gas is given by the 
Spitzer (1962) formula, 



1.84 X 10~T t 
In A 



5^^5/2 



(8) 



where Te is the electron temperature (Te = T for a one- 
temperature plasma) and In A is Coulomb logarithm that 
for T > 4.2 X W^K is 



lnA = 29.7 + Inn" ^/^(Te/lO^/f). 



(9) 



The heat is conducted by the electron, and equation (8) 
includes the effect of the self-consistent electric required to 
maintain the electric current at zero; this reduces k by a 
factor of about 0.4 from the value it would otherwise have 
(Cowie & McKee 1977, hereafter CM77). 

As noted in the introduction, the inner regions of hot 
accretion flows arc, in many cases, coUisionlcss with electron 
mean free path due to Coulomb collision larger than the ra- 
dius (e. g. Tanaka & Menou 2006). When the mean free 
path of an electron becomes comparable to or larger than 
the temperature gradient scale A > T/|VT|, equation (7) for 
the heat flux is no longer valid; CM77 described this effect 
as saturation. The maximum heat flux in a plasma can be 
expressed as (3 /2)nckT^Vchar , where Vchar is a characteristic 
velocity which one might expect to be the order of the elec- 
tron thermal velocity (Parker 1963). Assuming Maocwellian 
distribution for heat source, the characteristic velocity can 
be written as (Williams 1971; CM77) 



(10) 
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Similar to CM77, wc assume that the heat flux is reduced 
by the same factor of 0.4 in the saturated case as in the 
classical (coUisional) case so that the saturated heat flux is 



: OAn^kTe 
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CM77 showed that the saturated heat flux is signiflcantly 
less than conjectured by Parker (1963). Thus, in order to 
explicitly allow for uncertainty in the estimate of Fgat, they 
introduced a factor of 03, which was less than unity and 
rewrote equation (11) as 



■ 50s pCs 



(12) 



where Cs is sound speed, which is defined as = p/ p. The 
factor of (f)s is called as saturation constant (CM77). Now, 
the viscous heating rate and the energy transport by thermal 
conduction are expressed as 

2 



Qvis = vpr 
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dr 
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By using equations (13) and (14) for the advective transport 
of energy, we can write 

Q... = /.pr^(f)'-^|:(r^F.„.) (15) 

The mass accretion rate in a qausi-spherical accretion 
flow can be written as 



M{r,t) = — 47rr pVr 



(16) 



We will use this quantity in the next section, and will inves- 
tigate effects of saturation constant and viscous parameter 
on it. 



3 SELF-SIMILAR SOLUTIONS 
3.1 analysis 

Tanaka & Menou (2006) solved essentially equations (1)- 
(4) for the case of a steady, radially self-similar flow. Here, 
we will try to find unsteady self-similar solutions for these 
equations. Thus, we introduce a similarity variable f and 
assume that each physical quantity is given by the following 
form: 



i = r(GA^,^^)~^/^ 

p{r,t) = R{e,){Kh/GM,)t-\ 

p{r,t) = ^(e)(Mo/(GA/.)'/')^"^/^ 

«.(r,^) = ^/(C)(GM.)^/^^-l/^ 

n{r,t) = ujiOt'\ 
M{r,t) = MorhiO, 



(17) 
(18) 
(19) 
(20) 
(21) 
(22) 



where Mo is a constant and its value can be obtained by 
typical values of the system. In addition, wc assumed that 
M{r,t) under similarity transformations is a function of ^ 
only (Khesali & Faghei 2008, 2009). Substitution of above 
transformations into the basic equations (l)-(4), yields di- 
mensionless equations below. 
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These equations provide a fourth-order system of non-linear 
ordinary differential equations that must be solved numeri- 
cally. 



3.2 Inner limit 

An appropriate asymptotic solution as ^ — >■ is the form as 

j?(0~r'^'(ft) + i?iC + ---), (27) 

n(C)~r'^'(no + ni5 + ••■), (28) 

V{0-^r'^HVo + Vi^ + ---), (29) 
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(30) 



in which undetermined coefficients of Ro, Ri, Ho, and etc 
must be specified. By substituting above relations in equa- 
tions (23)-(26) and choosing the significant sentences, we 
can write 



Ho 

27?oV^o + 3ano « 0, 
6V^o(7-|) + (7-l) 



0, 



! 0. 



(31) 
(32) 

(33) 



Also, the dimensionless mass accretion rate, rri(^), under the 
above asymptotic solution becomes 



-AttRoVo, 



(34) 



where rhin is the value of rh at ^in, where is a point near 
to the centre. After algebraic manipulations for equations 
(31)- (34), we obtain an algebraic equation for Rq: 
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and the rest of the physical variables are 
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Without thermal conduction, (f>s = 0, equation (35) can be 
solved analytically. Since, we want to consider systems with 
non-zero saturation constant, 4>s 7^ 0, we will solve this equa- 
tion numerically. 

3.3 Numerical solution 

If the value of is guessed, i.e. we take a point very 
near to the centre, the equations (23)-(26) by Runge-Kuttar 
Fehlberg fourth-fifth order method can be integrated from 
this point outward by the above expansion [(27)- (30)]. Ex- 
amples of such solutions are presented in Figs 1-5. 

3.3.1 The influences of saturation constant and viscous 
parameter on physical quantities 

The delineated quantity of H/R in Figs 1 and 2 is the sound 
speed square in self-similar flow, which is rescaled in the 
course of time and represents the flow temperature. The pro- 
flies of n/_R in Fig. 1 show the flow temperature decreased 
by adding saturation constant, <j)s. Because, the generated 
heat by viscous dissipation can be transfered by thermal 
conduction. Furthermore, this temperature decrease is qual- 
itatively consistent with simulation results of Sharma et al. 
(2008) and Wu ct al. (2010). We know the viscous dissipa- 
tion of the flow increases by adding a parameter. Thus, the 
temperature increased by viscous parameter confirmed by 
the temperature profiles in Fig 2. The density profiles show 



the gas density increased by adding the (ps parameter. It 
can be due to temperature fall of fluid. Increasing density 
by adding saturation constant is another consistency of our 
results with simulations of Wu et al. (2010). Also, the den- 
sity proflles show that it decreases by adding the viscous 
parameter. This also could be result of the temperature ris- 
ing. The viscous turbulence in this paper is proportional 
to the gas temperature {v oc cl oc T). Thus, increase or 
decrease of temperature will affect dynamics of the accret- 
ing gas. As we said the temperature decreased by satura- 
tion constant implies that the viscous turbulence decreases, 
too. The decreasing of viscous turbulence reduces the effect 
of negative viscous torque in angular momentum equation. 
Thus, we expect the flow rotates faster by adding saturation 
constant that confirmed by the angular velocity profiles in 
Fig. 1. Since, the efficiency of the angular momentum trans- 
port decreases by adding the saturation constant, we expect 
the decrease of radial infall velocity that the radial veloc- 
ity profiles in Fig. 1 confirm it. The efficiency of angular 
momentum transport increases by adding the viscous pa- 
rameter of a. Thus, we expect the flow rotates slower and 
accretes faster by adding the viscous parameter. The profiles 
of radial and angular velocities in Fig. 2 confirm them. 



3.3.2 Mass accretion rate 

The behaviour of mass accretion rate as a function of simi- 
larity variable ^ for several values of the viscous parameter 
and saturation constant are plotted in Fig. 3. In the present 
model, the mass accretion rate is reduced by radius. While, 
the mass accretion rate in steady hot accretion flows is a con- 
stant (Tanaka & Menou 2006). There are some researches in 
steady hot accretion flows that have studied power-law func- 
tion of mass accretion rate (Shadmehri 2008; Abbassi et al. 
2008). However, the mass accretion rate in their models is 
not dependent on important parameters such as saturation 
constant and viscous parameter. The profiles of mass accre- 
tion rate in Fig. 3 show that it is reduced by adding the 
saturation constant. This property is qualitatively consis- 
tent with numerical results of Johnson & Quatacrt (2007). 
Also, the profiles of mass accretion rate imply that it in- 
creases by adding the viscous parameter of a. This property 
is qualitatively consistent with previous works in accretion 
flows (e. g. Park 2009). 



3.3.3 Mach number 

Here, it will be interesting to investigate the existence of the 
transonic point in hot accretion flow. The transonic point 
occurs in place that the amount of Mach number becomes 
equal to unity. The Mach number referring to the reference 
frame is defined as (Gaffet & Fukue 1983; Fukue 1984) 



Vr — Vf V — n$, 



where 

dr r 
Vf ~ —r = ri— 
dt t 



(39) 



(40) 



is the velocity of the reference frame which is moving out- 
ward as time goes by, which the sound speed can be subse- 
quently expressed as 
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Figure 2. Same as Figure 1, but </is = 0.03. The solid, dahsed, and short-dashed Hues represent a = 0.05, 0.1, and 0.2, respectively. 
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Figure 3. Time-dependent self-similar solution of mass accretion rate. The input parameters in left panel are same as Figure 1, but the 
solid, the dashed, and the short-dashed lines represent 4>s = 0.01, 0.02, and 0.03, respectively. The input parameters in right panel are 
same as Figure 2. but (ps = 0.01. 
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Figure 4. Time-dependent self-similar solution of Mach number. The input parameters in left panel are same as Figure 1, but a = 0.5 
and the solid, the dashed, and the short-dashed lines represent 4>s = 0.0, 0.2, and 0.3, respectively. The input parameters in right panel 
are same as Figure 1, but <j>s = 0.05 and the solid, the dashed, and the short-dashed lines represent a = 0.2, 0.3, 0.4, respectively. 



Cs = -^= S^{GM,/tf^^ (41) 

and, S = (n/R)^^^ the sound speed in self-similar flow is 
rescaled in the course of time. The Mach number introduced 
so far, represents the instantaneous and local Mach number 
of the unsteady self-similar flow. In steady self-similar solu- 
tion (e. g. Tanaka &: Mcnou 2006), the Mach number does 
not vary by radii and is a constant. While, the Mach num- 
ber in unsteady self-similar varies by radii (see Fig. 4). As 
seen in Fig. 4, there is a transonic point (|/i| = 1). The de- 
pendency of transonic point to saturation constant shows 
that this point moves inward by adding the parameter of 
(ps- Because, thermal conduction transfers the heat to larger 
radii, so the sound speed/temperature decreases by adding 
the saturation constant. In other words, whatever the radii 
smaller, the radial velocity relative to the sound speed larger. 
Also, the Mach number profiles show the transonic point de- 
creasing by adding the viscous parameter which can be due 
to increase of radial velocity along with adding the viscous 
parameter. 



3.3.4 Comparison of steady and unsteady self-sinular 
solutions 

As a comparison between steady and unsteady self-similar 
solution, the physical quantities in unsteady self-similar so- 
lutions are divided into their radial dependence in steady 
self-similar solution then plotted in terms of ^ in Fig. 5. Also, 
the effect of adiabatic index on hot accretion flow is inves- 
tigated in Fig. 5. The dchneatcd quantities {R^^^^, V 
■ ■ ■) iu Fig 5 arc constant in steady self-similar solutions of 
hot accretion flows (Tanaka & Menou 2006; Shadmehri 2008; 
Abbassi et al. 2008; Ghanbari et al. 2009). While they vary 
with position in this research. Fig. 5 represents the density 
profile varies shallower than ^ that this property is qual- 
itatively consistent with simulation results (e.g., Stone et al. 
1999; Igumcnshchev & Abramowicz 1999; Stone & Pringle 
2001; Hawlcy & Balbus 2002; Igumenshchev et al. 2003). The 
radial dependency of the temperature and radial velocity in 
unsteady self-similar solution show that they vary deeper 
than the steady self-similax solution. Also, the study of the 
angular velocity show that it varies shallower than 
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Figure 5. Same as Figure 1, but ■ 



: 0.1. The solid, the dashed, and the short-dashed lines represent 7 = 1.3, 1.4, and 1.5, respectively. 



Thus, the physical quantities in the present model avoid the 
limits of the steady self-similar solution. The physical quan- 
tities profiles in Fig. 5 show that their radial dependency in 
the unsteady self-similar limited to the steady self-similar 
solution by adding adiabatic index 7. 



4 SUMMARY AND DISCUSSION 

In hot accretion flows, the collision timescale between ions 
and electrons is longer than the inflow timescale. Thus, the 
inflow plasma is collisionless, and the transfer of energy by 
thermal conduction can be dynamically important. The low 
collisional rate of the gas is confirmed by direct observa- 
tion, particularly in the case of the Galactic centre (Quatacrt 
2004; Tanaka & Menou 2006) and in the intraclustcr medium 
of galaxy clusters (Sarazin 1986). 

Here, we have investigated how thermal conduction af- 
fects dynamics of hot quasi-spherical accretion flows. We 
adopted the presented solutions by Ogilvie (1999) and 
Tanaka & Menou (2006). Thus, we assumed that angular 
momentum transport is due to viscous turbulence and the 
a-prescription is used for the kinematic coefiicient of viscos- 
ity. We also assumed the flow does not have a good cooling 
efficiency and so a fraction of energy accretes along with 
matter on to the central object. The effect of thermal con- 
duction is studied by a saturation form of it introduced by 
Cowie & McKee (1977). To solve the equations that govern 



the dynamical behaviour of hot accretion flow, we have used 
unsteady self-similar solution. 

The effect of saturation constant and the viscous pa- 
rameter on the present model is investigated. The solutions 
show that with the increase of conductivity, the equatorial 
density becomes denser and the temperature becomes lower. 
These results are qualitatively consistent with simulation re- 
sults of Wu et al. (2010). Furthermore, the solutions show 
that by adding the saturation constant, the angular velocity 
becomes larger and the radial velocity decreases. The mass 
accretion rate is reduced by adding the saturation constant 
that is qualitatively consistent with the result of Johnson & 
Quataert (2007). The solutions imply that the viscous pa- 
rameter has opposite effects in comparison to saturation con- 
stant on physical quantities of the system. Also, the study 
of physical quantities of the present model in comparison 
to steady self-similar solution show that our results deviate 
from steady self-similar solution and do not have its limits. 

Here, we studied dynamical behaviour of hot accretion 
flow in one-dimensional approach ignored by latitudinal de- 
pendence of physical quantities. Although, some authors 
have shown that latitudinal dependence of physical quanti- 
ties is important for structure and dynamics of hot accretion 
flow (Tanaka & Menou 2006; Ghanbari et al. 2009; Wu et 
al. 2010). Thus, latitudinal behaviour of the present model 
can be investigated in other studies. 
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